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Weakly Interacting Massive Particles (WIMPs) are among the main candidates for 
the relic dark matter (DM). The idea of the direct DM detection relies on elastic spin- 
dependent (SD) and spin-independent (SI) interaction of WIMPs with target nuclei. 
In this review paper the relevant formulae for WIMP event rate calculations are 
collected. For estimations of the WIMP-proton and WIMP-neutron SD and SI cross 
sections the effective low-energy minimal supersymmetric standard model is used. 
The traditional one-coupling-dominance approach for evaluation of the exclusion 
curves is described. Further, the mixed spin-scalar coupling approach is discussed. 
It is demonstrated, taking the high-spin 73 Ge dark matter experiment HDMS as 
an example, how one can drastically improve the sensitivity of the exclusion curves 
within the mixed spin-scalar coupling approach, as well as due to a new procedure 
of background subtraction from the measured spectrum. A general discussion on the 
information obtained from exclusion curves is given. The necessity of clear WIMP 
direct detection signatures for a solution of the dark matter problem, is pointed out. 

PACS: 95.30.-k, 95.35.+d, 14.80.Ly, 12.60.Jv 

I. INTRODUCTION 

To our knowledge the galactic Dark Matter (DM) particles do not emit any detectable 
amounts of electromagnetic radiation and manifest themselves only gravitationally by affecting 
other astrophysical objects. The first evidence of this kind of substance came from the study 
of galactic rotation curves, i.e. from measurement of the velocity with which stars, globular 
stellar clusters, gas clouds, or dwarf galaxies orbit around their centers [1]. If the mass of 
these galaxies was concentrated in their visible parts, the orbital velocity at large radii r should 
decrease in accordance with Kepler's law as 1/y/r. Instead, it remains approximately constant 
to the largest radius where it can be measured. This implies that the total mass M(r) felt by 
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an object at a radius r must increase linearly with r. Studies of this type imply that 90% or 
more of the mass of the large galaxies is in their dark halos [2-4] . 

The mass density averaged over the entire Universe is usually expressed in units of the critical 
density p c ~ 10 _29 g/cm 3 . The dimensionless ratio Q = p/ p c = 1 corresponds to a flat Universe. 
Analyzes of galactic rotation curves imply Q > 0.1 (see for example, [5-7]). Studies of clusters 
and superclusters of galaxies through gravitational lensing or through measurements of their 
X-ray temperature, as well as studies of the large-scale streaming of galaxies favor larger values 
of the total mass density of the Universe f2 > 0.3 (see, for example [7, 8]). Finally, naturalness 
arguments and most inflationary models prefer fi = 1.0 to a high accuracy. The requirement 
that the Universe be at least 10 billion years old implies Qh 2 < 1, where h is the present 
Hubble parameter in units of 100 km/ (sec-Mpc) [7]. The total density of luminous matter only 
amounts to less than 0.4% of the critical density [9, 10]. Analyzes of Big Bang nucleosynthesis 
determine the total baryonic density to lie in the range 0.017 < il h h 2 < 0.024 [7]. The upper 
bound implies f2 b < 0.05, in obvious conflict with the lower bound f2 > 0.3. Most Dark Matter 
must therefore be non-baryonic. Some sort of "new physics" is required to describe this exotic 
matter, beyond the particles described by the Standard Model of particle physics. 

Exciting evidence for a flat and accelerating universe was claimed by [11-13]. The position 
of the first acoustic peak of the angular power spectrum (of the temperature anisotropy of the 
cosmic microwave background radiation [14, 15]) strongly suggests a flat universe with density 
parameter ft — 1 while the shape of the peak is consistent with the density perturbations 
predicted by models of inflation. Data support Q = Q M + Q A = 1 where Q M is the matter 
density in the universe and Q\ is usually assumed to be a contribution of a non-zero cosmological 
constant (the energy density of the vacuum). A first claim for the existence of a non- vanishing 
cosmological constant has been made already in 1986 [16, 17]. Recent investigations of the 
Cosmic microwave background temperature anisotropy by the Wilkinson Microwave Anisotropy 
Probe (WMAP) [18-20] and the galaxy power spectrum with the baryon acoustic peak by the 
Sloan Digital Sky Survey (SDSS) [21-23] supplied us with the values for the cosmological 
parameters given in Table I. The parameters unambiguously confirm the existence of a large 
amount of dark matter. We omit in this paper discussion of the Dark Energy — another 
mysterious substance which is connected with the accelerating Universe and fills [9] the gap 
between a flat Universe and the measured amount of Dark Matter (f2 DM + f2 DE = f2 tot = 1). 
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Hubble Constant 
Baryon Density 
Matter Density 


h = o.704±g:8ii 

n h h 2 = 0.0219 ±0.0007 
n M h 2 = 0.132 ±0.004 


Baryon/Critical Density 
Matter /Critical Density 
Total/Critical Density 
Age of the Universe 


n h = 0.0442 ± 0.003 
ft M = 0.249 ± 0.018 
n tot = 1.011 ± 0.012 
to = 13.7 ± 0.2 Gyr 



TABLE I: Some Basic Cosmological Parameters from WMAP and SDSS. From [7]. 



In 2006 an exiting "visualization" of the invisible dark matter substance (see Fig. 1) has been 
obtained by means of gravitational lensing [24] . 




FIG. 1: Three-dimensional reconstruction of the dark matter distribution. It is obtained by the 
Hubble Space Telescope Collaboration [24] from the differential growth of the gravitational lensing 
signal between many thin discrete redshift slices. The three axes correspond to right ascension, 
declination and redshift. The distance from the Earth increases toward the bottom of the picture. For 
details see the original paper [24]. 

According to the estimates, based on a detailed model of our Galaxy [25], the local density 
of DM (nearby the solar system) amounts to about p™ al ~ 0.3 GeV/cm 3 ~ 5 ■ 10~ 25 g/cm 3 , 
with an uncertainty within a factor of two [7]. It is assumed to have a Maxwellian velocity 
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distribution in the galactic rest frame with mean velocity v ~ 270 km/sec [26, 27]. The local 

flux of DM particles x is expected to be ~ — — ■ 10 5 cm _2 s _1 . This value is often 

m x 

considered as a promising basis for direct dark matter search experiments. 

Weakly Interacting Massive Particles (WIMPs) are among the most popular candidates for 
the relic dark matter. There is no room for such particles in the Standard model of parti- 
cle physics (SM). The lightest supersymmetric (SUSY) particle (LSP), the neutralino (being 
massive, neutral and stable) is currently often assumed to be a favorite WIMP dark matter 
particle. The nuclear recoil energy due to elastic WIMP-nucleus scattering is the quantity to 
be measured by a terrestrial detector in direct DM detection experiments [28]. Detection of 
the very rare events of such WIMP interactions is a challenge for modern particle physics, 
because of the very weak WIMP coupling with ordinary matter. The rates expected in SUSY 
models range from 10 to 10~ 7 events per kilogram detector material and day (see, for example 
[27, 29-38]). Moreover, for WIMP masses between a few GeV/c 2 and 1 TeV/c 2 , the energy 
deposited by the recoil nucleus is less than 100 keV. Therefore, in order to be able to detect 
a WIMP, an experiment with a low-energy threshold and an extremely low radioactive back- 
ground is required. Furthermore, to indeed detect a WIMP one has to unambiguously register 
some positive signature of WIMP-nucleus interactions (directional recoil or annual signal mod- 
ulation) [26, 27]. This means one has to perform a measurement with a detector of large target 
mass during several years under extremely low radioactive background conditions (see also the 
discussions of other complications in [39-41]). Despite of all these problems huge effort is at 
present put into direct detection of DM particles (see for example, [2, 7, 42, 43]). 

Till now only the DAMA (DArk MAtter) collaboration claims [44-47] observation of first 
evidence for a dark matter signal due to registration of the predicted annual modulation of 
specific shape and amplitude due to the combined motions of the Earth and Sun around the 
galactic center [26]. Aimed since more than one decade at the direct detection of DM particles, 
the DAMA experiment (DAMA/Nal) with 100 kg of highly radio-pure Nal(Tl) scintillator 
detectors successfully operated till July 2002 at the Gran Sasso National Laboratory (LNGS) 
of the I. N.F.N. On the basis of the results obtained over 7 annual cycles (107731 kg-day total 
DAMA exposure) the presence of a WIMP-model independent annual modulation signature 
was observed at a 6.3 a C.L. [46]. The main result of the DAMA observation of the annual 
modulation signature is a low-mass region of the WIMPs (40 < m x < 150 GeV) and relatively 
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high allowed SI or/and SD cross sections (for example, 1 • 1CT 7 pb < 0gj(O) < 3 ■ 1CT 5 pb), 
provided these WIMPs are cold dark matter particles. 

Although there are other experiments like EDELWEISS, CDMS, etc, which give sensitive 
exclusion curves, no one of them at present has the sensitivity to look for the modulation effect. 
Due to the small target masses and short running times these experiments are unable to see a 
positive annual modulation signature of the WIMP interactions. Some other experiments with 
much larger mass targets (mostly Nal) unfortunately are also unable to register the positive 
signature due to not good enough background conditions (see for example, [48-50]). Often the 
results of these and the DAMA experiment have been compared not on the basis of a complete 
analysis including simultaneously SI and SD WIMP nucleus interaction. This sometimes gives 
rise to quite some confusion in the literature (for a discussion see [51, 52]), and to attempts to 
reconcile an artificial DAMA "conflict" with the other experiments [53-59]. 

Despite of the well-known attempts of the DAMA collaboration to prove this observation 
with a new larger Nal setup DAMA/LIBRA [60], it is obvious that such a serious claim should be 
verified at least by one other completely independent experiment. To confirm this DAMA result 
one should perform a new experiment which would have (in reasonable time) the same or better 
sensitivity to the annual modulation signal (and also it would be reasonable to locate this new 
setup in another low-background underground laboratory). This mission, in particular, could 
be executed by new-generation experiments with large enough mass of germanium high purity 
(HP) detectors both with spin ( 73 Ge) and spin-less (natural Ge). Despite of obviously necessary 
strong figthing against backgrounds, the main direction in development of new-generation DM 
detectors concerns remarkable enlargement of the target mass to be able to observe these 
positive signatures, and thus to detect DM and to prove, or disprove the DAMA claim. In 
particular, an enlarged version of the EDELWEISS setup with 40 kg bolometric Ge detectors 
[61] together with, perhaps, SuperCDMS [62, 63], as well as enlarged ZEPLIN [64] or KIMS 
[65] experiments might become sensitive to the annual modulation in some future. 

The main efforts (and expectations) in present direct dark matter searches are concentrated 
in the field of so-called spin-independent (or scalar) interaction of a dark matter WIMP with 
a target nucleus. This is because it was found theoretically that for heavy enough nuclei 
this spin-independent (SI) interaction of DM particles with nuclei usually gives the dominant 
contribution to the expected event rate of its detection. The reason is the strong (proportional 
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to the squared mass of the target nucleus) enhancement of the SI WIMP-nucleus interaction. 

The spin-1/2 WIMP particles, like the LSP neutralinos, interact with ordinary matter pre- 
dominantly by means of axial vector (spin-dependent) and vector (spin-independent) couplings. 
There is some revival of interest in the WIMP-nucleus spin-dependent interaction from both 
theoretical (see e.g. [35, 36, 52, 66-70]) and experimental (see e.g. [57, 71-79]) points of view. 
There are some proposals aimed at direct DM detection with relatively low-mass isotope tar- 
gets [71, 72, 76-78, 80] as well as some attempts to design and construct a DM detector which 
is sensitive to the nuclear recoil direction [81-87]. Low-mass targets make preference for the 
low-mass WIMPs and are more sensitive to the spin-dependent WIMP-nucleus interaction as 
well [29, 35, 66, 68, 70, 88, 89]. 

There are at least three reasons to think that SD (or axial-vector) interaction of the DM 
WIMPs with nuclei could be very important. First, contrary to the only one constraint for SUSY 
models available from the scalar WIMP-nucleus interaction, the spin WIMP-nucleus interaction 
supplies us with two such constraints (see for example [68] and formulae below). Second, one 
can notice [35, 36] that even with a very sensitive DM detector (say, with a sensitivity of 
10~ 5 events/day/kg) which is sensitive only to the WIMP-nucleus scalar interaction (with spin- 
less target nuclei) one can, in principle, miss a DM signal. To safely avoid such a situation 
one should have a spin-sensitive DM detector, i.e. a detector with spin-non-zero target nuclei. 
Finally, there is a complicated nuclear spin structure, which possesses the so-called long g-tail 
form-factor behavior. The SI WIMP-nucleus cross section, despite being proportional to A 2 , 
vanishes very quickly (exponentially) with increasing momentum transfer q 2 . The SD WIMP- 
nucleus cross section decreases not so quickly with q 2 and remains still final at the recoil energies 
(En = q 2 /(2M4)), where the SI cross section is already zero. Therefore for heavy mass target 
nuclei and heavy WIMP masses the SD efficiency to detect a DM signal could be much higher 
than the SI efficiency [66]. Therefore, simultaneous study of both spin-dependent and spin- 
independent interactions of the DM particles with nuclei significantly increases the chance to 
observe the DM signal [36, 51, 69, 90]. 

Following R.Bernabei et al. [45, 46] it was stressed in [51, 52] that for analyzing data 
from DM detectors with spin-non-zero targets one should use the so-called mixed spin-scalar 
coupling approach. This approach is used to demonstrate, taking the high-spin 73 Ge detector 
HDMS [91, 92] as example, how one can stronger improve the exclusion curves. The mixed 
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spin-scalar coupling approach allowed one to extract information about both SI and SD WIMP- 
nucleon cross sections analyzing background spectra from the two HDMS setups (prototype and 
final) simultaneously. This procedure allows an improvement (see our new analysis in [93]) of 
the exclusion curves relative to the relevant curves obtained in the traditional one-coupling 
dominance approach for the HDMS in [79]. 

The present paper is organized as follows. In the next Section the main formulae for event 
rate calculations are collected. In Section III the effective low-energy minimal supersymmetric 
standard model (effMSSM) is used for calculation of the WIMP-proton and WIMP-neutron 
SD and SI cross sections. In Section IV the traditional one coupling dominance approach for 
evaluation of the exclusion curves is discussed. In Section V the mixed spin-scalar couplings 
approach is described, the DAMA-inspired exclusion domains for both above-mentioned cou- 
plings are given and compared with SUSY calculations. In Section VI the mixed spin-scalar 
coupling scheme is applied to the high-spin 73 Ge dark matter search experiment HDMS. It is 
demonstated how one can strongly improve the quality of the exclusion curves within the mixed 
spin-scalar coupling approach as well as by using a new procedure of background subtraction 
from the measured spectrum. In Section VII a general discussion is given. The conclusion 
summarizes the main items of this review paper. 



Many experiments try to detect directly a relic DM WIMP (or neutralino) x with mass m x 
via its elastic scattering on a target nucleus (A, Z). The nuclear recoil energy E R is measured 
by a proper detector deeply underground (Fig. 2). The differential event rate in respect to 
the recoil energy (the spectrum) is the subject of the measurements. The rate depends on the 
density and the velocity distribution of the relic WIMPs in the solar vicinity f(v) and the cross 
section of WIMP-nucleus elastic scattering [27, 29, 33, 38, 89, 94-96]. The differential event 
rate per unit mass of the target material has the form 



We assume these WIMPs (or neutralinos) to be the dominant component of the DM halo of our 
Galaxy with a density p x = 0.3 GeV/cm 3 in the solar vicinity. The (real) nuclear recoil energy 
-E'r — <? 2 /(2M4) is typically about 10 _6 m x and Nt is the number density of target nuclei with 



II. EVENT RATE AND CROSS SECTIONS 
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FIG. 2: Detection of the cold dark matter (WIMPs) by elastic scattering from target nuclei in the 
detector. Due to the expected annual modulation signature of the event rate Eq. (1) the Sun-Earth 
system is a particularly proper setup for successful direct DM detection. 

1 /2 

mass Ma- v mSLX = v esc ~ 600 km/s, v min = (MaEr/2/j<\) is the minimal WIMP velocity 
which still can produce the recoil energy Er. The WIMP-nucleus differential elastic scattering 
cross section for spin-non-zero (J ^ 0) nuclei contains coherent (spin- independent, or SI) and 
axial (spin-dependent, or SD) terms [66, 97, 98]: 

The normalized (-^sdsiW = -0 finite- momentum-transfer nuclear form-factors -^sdsi^ 2 ) = 
— ^ — -— can be expressed through the nuclear structure functions as follows [66, 97, 98]: 

S&(q) = E \(J\\C L (q)\\J)\ 2 ^\(J\\C ( q )\\J}\ 2 , 

L even 

SUt) = E (\(N\\T£ l5 (q)\\N}\ 2 + \(N\\Cl(q)\\N}\ 2 ). (3) 

Lodd 

The explicit form of the transverse electric T el5 (q) and longitudinal C 5 (q) multipole projections 
of the axial vector current operator and the scalar function Cl(?) can be found in [66, 70, 97-99]. 
For q = the nuclear SD and SI cross sections can be presented as follows 
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^(0) = f^f ^^K(S^) + a„(S-)} 2 (5) 

^D(0){(S^)co S ^+(S^)sin^} 2 . (6) 



sm "' (2J 

/4 4 J + 1 



fj%3 J 

Following Bernabei et al. [45, 46] the effective spin WIMP-nucleon cross section o"sd(0) and 
the coupling mixing angle 6 were introduced 



a,. 



= ^< + <, tan£ = ^; (7) 



7T 3 L " " J dp 



^SD = ^SD • COS ^SD = ^SD • sin 6 - ( 8 ) 



Here, a a = — ;rr- is the reduced mass of the neutralino and the nucleus, and it is assumed 

' ^ A m x + M A 

that f/,^ = fi 2 . The dependence on effective WIMP-quark (in SUSY neutralino-quark) couplings 
C q and A q in the underlying theory 

£eff = Yl (A - X7 M 75X • 97^759 + C q ■ XX • 5?) + ••• (9) 

and on the spin (A<f and the mass or scalar (ffl structure of the proton and neutron 

enter into these formulae via the zero-momentum-transfer WIMP-proton and WIMP-neutron 
SI and SD cross sections: 

,2 



<(0) = 4^4 c = eg'" = £CJ(^); (10) 

7T 

o^(o) = i24x„ « P = EM P) , «n = EAA< n) . (ii) 



The factors A q p ' n \ which parameterize the quark spin content of the nucleon, are defined as 
2/\( n 'P) s ^ = (j,^ slipq^^ipqlp, s)( p ,n). The quantity (S^^) denotes the total spin of protons 
(neutrons) averaged over all A nucleons of the nucleus (A, Z): 

(S^ B) > = (A\tf {n) \A) = (A\J2s; {n) \A). (12) 

i 

The mean velocity (v) of the relic neutralinos of our Galaxy is about 300 km/s = 10~ 3 c. 
Assuming q max R <C 1, where R is the nuclear radius and g max = 2/i^f is the maximum of the 
momentum transfer in the process of the xA scattering, the spin-dependent matrix element 
takes a simple form (zero momentum transfer limit) [100, 101]: 

M = C{A\a p S p + a n S n \A) ■ s x = CA(A\3\A) • s x . (13) 
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Here, s x denotes the spin of the neutralino, and 

_ (N\a p S p + a n S n \N) _ (N\(a p S p + a n S n ) ■ J\N) 

(N\J\N) ~ J(J+1) _ J J ' 1 J 

Note a coupling of the spin of x to the spin carried by the protons and the neutrons. The 
uncertainties arising from the electroweak and QCD scale physics are incorporated in the factors 
a p and a n . The normalization factor C involves the coupling constants, the masses of the 
exchanged bosons and the mixing parameters relevant to the lightest supersymmetric particle 
(LSP), i.e., it is not related to the associated nuclear matrix elements [102]. In the limit of zero 
momentum transfer q = the spin structure function in Eq. (3) reduces to the form 



S A (0) = h(A\\Y: ka + a.T^WA)? = ^±±J{J + 1)A 

47T „■ Z 7T 
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For the most interesting isotopes either (S^) or (S^) dominates ((S^) <C (S^)). See, for 
example, Table II. 

The differential event rate (1) can be given also in the form [46, 51]: 

dR(E R ) 



dE R 



k si (E r , m x ) cr S i + %d(^R> m x) °sd- (15) 



Ksi(^,™x) = ^|^3si(£r) 



K SD (E R ,m x ) = N T ^p^B SD (E R ) l^±l((S p )cos^+(S„)sin^) 2 



2m x n 



(16) 



B S i,sb(E r ) = -j-^rFg 1SD (E R )I(E R ). 



(v 2 

The dimensionless integral I(E R ) is a dark-matter-particle velocity distribution correction (see 
Eq. (15)): 

I(Er) = { ~fl I tMdx = ^ 3 + 2 f [erf^+^-erf^-??)], (17) 

(v) Jx min v 2 y/ir(l + 2ri 2 )ert(ri) + 2r]e " 

where one assumes that in the rest frame of our Galaxy WIMPs have a Maxwell-Boltzmann 

velocity distribution, and uses the dimensionless Earth speed with respect to the halo rj, as well 

3 Ad E 2 f x 2 

as x 2 ■ = [26, 27]. The error function is erf(x) = —= / dte~ l . The velocity variable 

4 [i\v z y 7r Jo 

is the dispersion v ~ 270km/c. The mean WIMP velocity (v) = \J\v. We also assume both 
form-factors F§ lSD (E R ) in the simplest Gaussian form following [103, 104]. In particular, this 
allows rather simple formulae (see Eq. (16)) to be used. Integrating the differential rate Eq. (1) 
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TABLE II: Zero momentum spin structure of nuclei in different models. The measured magnetic 
moments used as input are enclosed in parentheses. The variation of the (Sp ) and (S^) for fixed A 
reflects the level of inaccuracy and complexity of the current nuclear structure calculations. From [70]. 



M Ge (Lj = G 9/2 ) 


(S P ) 


ICS \ 

(Sn) 


fi (in hn) 


ISPSM, Ellis-Flores [103, 104J 





0.5 


— 1.913 


OGM, Engel-Vogel [105] 





0.23 


/ t~~\ O \ 

(-0.879) cxp 


1B1M, lacheilo et al. [106] and [98] 


—0.009 


0.469 


— 1.785 


IBFM (quenched), lacheilo et al. [106] and [98] 


-0.005 


0.245 


(-0.879) cxp 


TFFS, Nikolaev-Klapdor-Klemgrothaus, [107] 





0.34 




SM (small), Ressell et al. [98] 


0.005 


0.496 


-1.468 


SM (large), Ressell et al. [98] 


0.011 


0.468 


-1.239 


Oiv i\~ / 1 1 1\ i ~» n j_ i l"r\ol 

SM (large, quenched), Ressell et al. [98] 


0.009 


0.372 


(-0.879) cxp 


Hybrid bM, Uimitrov et al. [lOoJ 


U.UoU 


0.378 


-0.920 


127t It 7-1 A 

I (£j = #5/2) 


( s p) 


(Sn) 


fi (in fi N ) 


1SFSM, Ellis-1 lores [104, 109] 


1/2 





4.793 


OGM, Engel-Vogel [105] 


0.07 





(2.813) cxp 


ID T \ "\ T T 1 11 i 1 r -1 n^l 

IBFM, lacheilo et al. [106] 


r\ AHA 

0.464 


0.010 


(2.813) ex p 


T"n"Pl\ J / 1 1\ T 1 11 J_ 1 l~1 A/"l 

IBFM (quenched), lacheilo et al. [106] 


0.154 


0.003 


(2.813) ex p 


TFFS, Nikolaev-Klapdor-Klemgrothaus, [107] 


Air 

0.15 







SM (Bonn A), Ressell-Dean [101] 


0.309 


0.075 


2.775 {2.470} cff 


SM (Nijmegen II), Ressell-Dean [101] 


0.354 


0.064 3.150 {2.7930} eff 


131 Xe (Lj = £> 3/2 ) 


(Sp) 


(Sn) 


Ii (in n N ) 


ISPSM, Ellis-Flores [103, 104] 





-0.3 


1.148 


OGM, Engel-Vogel [105] 


0.0 


-0.18 


(0.692) cxp 


IBFM, lacheilo et al. [106] 


0.000 


-0.280 


(0.692) cxp 


IBFM (quenched), lacheilo et al. [106] 


0.000 


-0.168 


(0.692) cxp 


TFFS, Nikolaev-Klapdor-Kleingrothaus, [107] 




-0.186 




SM (Bonn A), Ressell-Dean [101] 


-0.009 


-0.227 


0.980 {0.637} cff 


SM (Nijmegen II), Ressell-Dean [101] 


-0.012 


-0.217 


0.979 {0.347} cfr 


QTDA, Engel [66] 


-0.041 


-0.236 


0.70 



from the recoil energy threshold e to some maximal energy e one obtains the total detection 
rate R(e, e) as a sum of the SD and SI terms: 

R(e,e) = R S i(e,e) + R SB (e,e) = dE R K sl (E R ,m x ) cr SI + J dE R K SD (E R ,m x ) a SB . (18) 

To accurately estimate the event rate i?(e, e) one needs to know a number of quite uncertain 
astrophysical and nuclear structure parameters as well as the very specific characteristics of an 
experimental setup (see, for example, discussions in [46, 110]). 



12 



As m x increases, the product qR starts to become non-negligible and the finite momentum 
transfer limit must be considered [70, 97-99, 101]. The formalism is a straightforward extension 
of that developed for the study of weak and electromagnetic semi-leptonic interactions in nuclei 
[98, 101]. With the isoscalar spin coupling constant = a n + a p and the isovector spin coupling 
constant a\ = a p — a n one can split the nuclear structure function S (q) (from Eqs. (2) and 
(3)) into a pure isoscalar term, S^(q), a pure isovector term, S^(q), and an interference term, 
S^q), in the following way: 

S A (q) = a 2 X (q) + ajS^q) + a a 1 S A 1 (q). (19) 

The relations S A (0) = C(J)«S P ) + (S n » 2 , S£(0) = C(J)«S P ) - (S n » 2 , and S&(0) = 

2C(J)((S 2 ) — (S 2 )) with C(J) = — ^— , connect the nuclear spin structure function 

S A (q = 0) with proton (S p ) and neutron (S n ) spin contributions averaged over the nucleus [99]. 

These three partial structure functions S^(q) allow calculation of spin-dependent cross sec- 
tions for any heavy Majorana particle as well as for the neutralino with arbitrary composition 
[100]. 

The first model to estimate the spin content in the nucleus for the dark matter search was the 
independent single-particle shell model (ISPSM) used originally by Goodman and Witten [28] 
and later in [94, 103, 111]. There are several approaches to more accurate calculations of the 
nuclear structure effects relevant to the dark matter detection. The list of the models includes 
the Odd Group Model (OGM) of Engel and Vogel [105] and their extended OGM (EOGM) 
[97, 105]; Interacting Boson Fermion Model (IBFM) of Iachello, Krauss, and Maino [106]; 
Theory of Finite Fermi Systems (TFFS) of Nikolaev and Klapdor-Kleingrothaus [107]; Quasi 
Tamm-Dancoff Approximation (QTDA) of Engel [66]; different shell model treatments (SM) 
by Pacheco and Strottman [112]; by Engel, Pittel, Ormand and Vogel [113] and Engel, Ressell, 
Towner and Ormand, [100], by Ressell et al. [98] and Ressell and Dean [101]; by Kosmas, 
Vergados et al. [31, 88, 114]; the so-called "hybrid" model of Dimitrov, Engel and Pittel [108] 
and perturbation theory based on calculations of Engel et al. [100]. For the experimentally 
interesting nuclear systems 29 Si and 73 Ge very elaborate calculations have been performed by 
Ressell et al. [98]. In the case of 73 Ge a further improved calculation by Dimitrov, Engel 
and Pittel was carried out [108] by suitably mixing variationally determined triaxial Slater 
determinants. At the present time the necessity for more detailed calculations especially for 
the spin-dependent component of the cross sections for heavy nuclei is well motivated. 
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To perform modern data analysis in the finite momentum transfer approximation it looks 
reasonable to use formulae for the differential event rate (Eq. (1)) as schematically given below: 



Here the ratio of isovector-to-isoscalar nucleon couplings is u = ai/a . The detector threshold 
recoil energy e and the maximal available recoil energy e (e < E R < e) have been introduced 
already in Eq. (18). In practice, for example with an ionization or scintillation signal, one has 
to take into account the quenching of the recoil energy, when the visible recoil energy is smaller 
than the real recoil energy transmitted by the WIMP to the target nucleus. 

Formulae (20) allow experimental recoil spectra to be directly described in terms of only 
three [68] (it is rather reasonable to assume Cgj(0) ~ °"si(0)) independent parameters (ogj, 
dp and a;) for any fixed WIMP mass m x and any neutralino composition. Comparing this 
formula with the observed recoil spectra for different targets (Ge, Xe, F, Nal, etc) one can 
directly and simultaneously restrict both isoscalar c (via <7gj) and isovector neutralino-nucleon 
effective couplings ao,i- These constraints, based on the nuclear spin structure functions for 
finite q, will impose the most model-independent and most accurate restrictions on any SUSY 
parameter space. Contrary to some other possibilities (see, for example, [46] and [115]), this 
procedure is direct and uses as much as possible the results of the accurate nuclear spin structure 
calculations. 

It is seen from Eqs. (8) and (20) that the SD cross sections <7g D and <7g D (or equivalently a\ 
and oj = ai/ao) are the only two WIMP-nucleon spin variables which can be constrained (or 
extracted) from DM measurements. Therefore there is no sense to extract from the data (with 
"artificial" twofold ambiguity) effective WIMP-nucleon couplings a p and a n . 

III. CROSS SECTIONS IN THE EFFECTIVE LOW-ENERGY MSSM 

To estimate the expected direct DM detection rates (with formulae (1), (18) or (20)) one 
should calculate cross sections <7 SI and a SD (or WIMP-nucleon couplings c and a p n ) within 



dR(e, e) 
dE R 

J\f(e,e,E R ,m x ) 



W(e, e, E R , m x ) [v S i(Er, m x ) + r]' SB (E R , m x , u) ajj] ; 




(20) 



r]' SB (E R ,m x ,u;) 



r} SI (E R ,m x ) 



{^ s 2 i(£r)}; 
^{w+i ^ oo(g) + u2 Su{q) + u Soi(q) )}- 
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the framework of some SUSY-based theory or take them from some experimental data (if it is 
possible). 

To obtain as much as general SUSY predictions it appeared more convenient to work within 
a phenomenological effective low-energy minimal SUSY model (effMSSM) whose parameters 
are defined directly at the electroweak scale, relaxing completely constraints following from any 
unification assumption (see for example, [8, 30, 33, 35, 36, 38, 69, 89, 116-124]). The effMSSM 
parameter space is determined by entries of the mass matrices of neutralinos, charginos, Higgs 
bosons, sleptons and squarks. The list of free parameters includes tan/5, the ratio of neutral 
Higgs boson vacuum expectation values; /i, the bilinear Higgs parameter of the superpotential; 
M li2 , soft gaugino masses; Ma, the CP-odd Higgs mass; m~, m~, m 2 ~ (m~, m'~), squark 
(slepton) mass parameters squared for the 1st and 2nd generation; m~ , rri~, m 2 ~ (m~ , m~), 
squark (slepton) mass parameters squared for the 3rd generation; A t , A^, A T , soft trilinear 
couplings for the 3rd generation. The third gaugino mass parameter M 3 defines the mass of 
the gluino in the model and is determined by means of the GUT assumption M 2 = 0.3 M 3 . 
In the MSSM the lightest neutralino x = Xi i s a mixture of four superpartners of gauge 
and Higgs bosons (Bino, Wino and two Higgsinos): x — NnB° + N 12 W° + N 13 H® + N^H®. 
Concerning the neutralino composition it is commonly accepted that x is mostly gaugino-like 
if P = + N 2 2 > 0.9 and Higgsino-like if P < 0.1, or mixed otherwise. 

To constrain the huge effMSSM parameter space and to have reliable predictions for the 
dark matter experiments one usually takes into account available information from colliders, 
astrophysics and rare decays. In our previous considerations [35, 36, 69, 121-124] the exper- 
imental upper limits on sparticle and Higgs masses from their non-observations [7, 125] were 
included. Also the limits on the rare b — > decay [126, 127] following [128-131] have been 
imposed. 

Furthermore, for each point in the MSSM parameter space (MSSM model) the relic density of 
the light neutralinos Vl x h 2 was evaluated with the code [121-123] based on the code DarkSUSY 
[132] with the allowance for all coannihilation channels with two-body final states that can 
occur between neutralinos, charginos, sleptons, stops and sbottoms as long as their masses are 
rrii < 2m x . Two cosmologically interesting regions were considered. One is 0.1 < Q x h 2 < 0.3 
and the other is the WMAP-inspired region 0.094 < Vt x h 2 < 0.129 [18, 19]. The possibility the 
LSP to be not the only DM candidate, with much smaller relic density 0.002 < Vlh 2 < 0.1 is 
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also taken into account. Further details can be found in [51]. In numerical studies of [35, 121- 
123, 133] the parameters of the effMSSM are randomly varied in the following intervals: 



-1 TeV < Mi < 1 TeV, -2 TeV < M 2 , n,A t <2 TeV, 
l<tan/3< 50, 60 GeV < Ma < 1000 GeV, 
10 GeV 2 < m 2 Q ,m 2 L ,m 2 Q3 ,m 2 L3 < 10 6 GeV 2 . 



.21) 



For the other sfermion mass parameters as before in [35, 36, 69, 121-124] we used the relations 



mi 



m|, m| 



2 2 

m~, m~ 
V T 



m~ = m 2 ~ , m 2 ~ = m~ . The parameters Ah and A T are 
Q 3 ' e 3 l 3 F 6 



c/ d Q /; /. / /; 
fixed to be zero. We consider the domain of the MSSM parameter space, in which we perform 

our scans, as quite spread and natural. 

Typical WIMP-nucleon cross sections of both spin (SD) and scalar (SI) interactions as 

function of the WIMP mass are shown as scatter plots in Fig. 3. In the figure open circles 
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FIG. 3: Cross sections of spin-dependent and spin-independent interactions of WIMPs with proton 
and neutron. Filled triangles (open circles) correspond to relic neutralino density 0.1 < ^ x /i 2 ^ 0.3 
(0.025 < Q x h 2 < 1). From [35, 121-123, 133]. 
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correspond to cross sections calculated under the old assumption that 0.025 < Q, x h 2 < 1. Filled 
triangles give the same cross section but the constraint on the flat and accelerating universe 
is imposed by 0.1 < VL x h 2 < 0.3. One can see that the reduction of the allowed domain for 
the relic density does not significantly affect spin-dependent and the spin-independent WIMP- 
nucleon cross sections. The different behavior of SD and SI cross sections with mass of the LSP 
can be seen from the plots. There is a more stringent lower bound for the spin-dependent cross 
section. It is at a level of 10" 7 pb. 

For more accurate investigation of the DAMA-inspired domain of the lower masses of the 
LSP (m x < 200 GeV) in [51] both a SD and <r SI have also been calculated within the effMSSM. 
To this end the intervals of the randomly scanned MSSM parameter space in [51] were narrowed: 

-200 GeV <M l < 200 GeV, -1 TeV < M 2 , n < 1 TeV, 50 GeV < M A < 500 GeV. (22) 

The results of these evaluations are shown as scatter plots in Fig. 4, which is the WIMP low- 
mass update of Fig. 3. In Fig. 4 filled green circles correspond to cross sections calculated when 
the neutralino relic density should just not overdose the Universe (0.0 < ^l x hl < 1.0). Filled red 
squares show the same cross sections when one assumes the relic neutralinos to be not the only 
DM particles and give only a sub-dominant contribution to the relic density 0.002 < J1 x /iq < 0.1. 
In the left panel of Fig. 4 these cross sections are shown with the black triangles corresponding 
to the case when the relic neutralino density is in the bounds previously associated with a 
so-called flat and accelerating Universe 0.1 < Q x hl < 0.3. The black triangles in Fig. 4 (right 
panel) correspond to the the WMAP and SDSS [18, 19] constraint on the matter relic density 
0.094 < £l x hl < 0.129 imposed in 2004. Despite a visible reduction of the allowed domain for 
the relic density due to the WMAP+SDSS result the upper bounds for the spin-dependent and 
the spin-independent WIMP-nucleon cross section are not significantly affected. 

Finally it is perhaps the right place here, to comment the following. Unfortunately the 
MSSM parameter space is huge and to obtain some reliable feeling, concerning, for example, 
the expected rate of dark matter detection when all relevant experimental and cosmological 
constraints are taken into account, one has nothing but this statistical numerical method (see 
for example, [29, 38, 89, 95, 96, 119, 134, 135]). This method allows lower and upper bounds 
for any observable to be estimated, and to make conclusions about the prospects for dark 
matter detection with present or future high-accuracy dark matter detectors. The larger the 
amount of points which confirm such a conclusion the better. The conclusions we made here are 
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FIG. 4: Cross sections of the spin-dependent (spin) and the spin-independent (scalar) interactions 
of WIMPs with the proton and the neutron. Filled green circles correspond to the relic neutralino 
density < ^x^o < 1> re d squares correspond to the sub-dominant relic neutralino contribution 
0.002 < O x /tg < 0.1 and black triangles correspond to the relic neutralino density 0.1 < ^yh^ < 0.3 
(left panel) and to the WMAP relic density 0.094 < n x h^ < 0.129 (right panels). 

based on hundreds of thousands of points which passed all constraints. Of course, we have no 
proved protection against peculiar choices of parameters which could lead to some cancellation 
and to small cross sections even if Higgs masses are small. Nevertheless, the probability of 
these choices is very small (about 1/100000), otherwise we should already meet them with our 
random scanning. On the other side, if these peculiar choices exist and one day would manifest 
themselves, this would be a very interesting puzzle, because it would be some kind of fine tuning 
of parameters, which requires strong further development of our understanding of the theory 
[124]. 
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IV. ONE-COUPLING DOMINANCE APPROACH 

From the definitions of SD and SI WIMP-nucleus and WIMP-nucleon cross sections 
(Eqs. (4)-(8), (10) and (11)) one can conclude that the spin observables in DM search give 
us two independent constraints on a SUSY model via 0g D (O) and o"g D (0), or, equivalently, via 
a p and a n (or a and a-J. These constraints are usually presented in the form of exclusion 
curves obtained with different target nuclei and recalculated in terms of nuclear-independent 
<7g D (0) (see for example, Fig. 5) and crg D (0) (see for example, Fig. 6). For a fixed mass of the 
WIMP the cross sections of SI or SD elastic WIMP-nucleon interaction located above these 
curves are excluded. 



D. 




FIG. 5: Exclusion curves for the spin-dependent WIMP-proton cross sections function of 

the WIMP mass. All curves, except the NAIAD and Tokio-LiF, are obtained in the one-coupling 
dominance approach with <r SI = and o"g D = 0. The DAMA/NaI-a(f) contours for the WIMP-proton 
SD interaction in 127 I are obtained on the basis of the positive signature of annual modulation within 
the framework of the mixed scalar-spin coupling approach [45, 46]. For details see [51]. 
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FIG. 6: Exclusion curves for the spin-dependent WIMP-neutron cross sections (o"g D versus the WIMP 
mass). The DAMA/NaI-a(f) contours for the WIMP-neutron SD interaction (subdominating in 127 I) 
are obtained from the relevant figures of [45, 46]. Note that the NAIAD curve [136] here corresponds to 
the WIMP-neutron SD interaction subdominant for 127 I. The WIMP-proton SD interaction dominates 
for this nucleus. The curve was obtained in the approach of [115]. It is much weaker in comparison 
with both the DAMA/Xe and HDMS-2003 curves. (For more details see [51] and Fig. 12). 

This simple presentation allows one to compare directly sensitivities of DM experiments with 
different nuclear targets. At the current level of accuracy (when « n n ) and 0^(0) ~ <Jgj(0), 
see Fig. 3) there is only one constraint for a WIMP-nucleon cross section (see Fig. 7) from 
spin-independent DM search experiments. Indeed, for the spin-zero nuclear target the exper- 
imentally measured event rate (Eq. (1)) of direct DM particle detection, via formula (2) is 
connected with the zero-momentum WIMP-proton (neutron) cross section (4). The zero mo- 
mentum scalar WIMP-proton(neutron) cross section crgj(O) can be expressed through effective 
neutralino-quark couplings C q (9) by means of expression (10). These couplings C q (as well as 
A q ) can be directly connected with the fundamental parameters of a SUSY model such as tan/3, 
Mi )2 , n, masses of sfermions and Higgs bosons, etc. Therefore experimental limitations on the 
SI neutralino-nucleon cross section supply us with a constraint on the fundamental parameters 
of an underlying SUSY model. In the case of the SD WIMP-nucleus interaction from a mea- 
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sured differential rate Eq. (1) one first extracts a limitation for o"g l D (0), and therefore has in 
principle two constraints [68] for the neutralino-proton a and neutralino-neutron a n effective 
spin couplings, as follows from relation (5). From Eq. (5) one can also see that contrary to 
the SI case (4) there is, in general, no factorization of the nuclear structure for cr^ D (0). Both 
proton (Sp ) and neutron (S^) spin contributions simultaneously enter into formula (5) for the 
SD WIMP-nucleus cross section o"^ D (0). Nevertheless, for the most interesting isotopes either 
<S£> or (S£) dominates ((S^) « (S£ b) )) [51, 70]. 




FIG. 7: WIMP-nucleon cross section (ogj(0)) limits in pb for spin-independent (scalar) interactions as 
a function of the WIMP mass in GeV. Shown are contour lines for some of the present experimental 
limits (solid lines) and some of projected experiments (dashed lines). All curves are obtained in 
the one-coupling dominance approach with <7 SD = 0. For example, the closed DAMA/Nal contour 
corresponds to complete neglection of SD WIMP-nucleon interaction. The open DAMA contour is 
obtained in [46] with the assumption that <7 SD = 0.08 pb > 0. Theoretical expectations for Ogj in the 
effMSSM from [51] are also shown by scatter plots for a relic neutralino density 0.09 < Q, x hQ < 0.13 
(black triangles) (see [51]). 

In earlier considerations [28, 94, 97, 103, 105, 111] one reasonably assumed that the nuclear 
spin was carried by the "odd" unpaired group of protons or neutrons and only one of either (S^) 
or (Sp ) was non-zero. In this case all possible non-zero-spin target nuclei can be classified into 
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n-odd and p-odd groups. Following this classification, the current experimental situation for 
the spin-dependent WIMP-proton and WIMP-neutron cross sections is naturally presented 
separately in Fig. 5 and Fig. 6. The DAMA/NaI-a(f) contours for the WIMP-proton SD 
interaction (dominating in 127 I) obtained on the basis of the positive signature of the annual 
modulation (closed contour) [46] and within the mixed coupling framework (open contour) 
[45] are also presented in Fig. 5. Similarly, the DAMA/NaI-a(f) [46] contours for the WIMP- 
neutron SD interaction (subdominant in 127 I) are given in Fig. 6. There are also exclusion 
curves for the SD cross section from the CDMS [137] and EDELWEISS [138] experiments with 
natural-germanium bolometric detectors (due to the small Ge-73 admixture). 

To compare experimental data with theoretical estimations in the effMSSM [51] one can 
superimpose the scatter plots for the SD and SI LSP-proton and LSP-neutron cross sections 
(from Fig. 4 or Fig. 3) in Figs. 5, 6, and 7. This is a traditional way to perform this comparison. 
One can easily see that both calclulated SD LSP-proton and LSP-neutron cross sections fall 
below the frames of Figs. 5 and 6, respectively. In particular, this means that experimental 
data (available at present in the form of exclusion curves) do not allow one to restrict the 
SUSY LSP-nucleon spin couplings. This is not the case for the SI WIMP-nucleon coupling. 
The scattered points (black triangles) for Ugj calculated in the effMSSM are clear seen in Fig. 7. 
Some of these points are already excluded by the DM measuremnts. 

Nevertheless, one would like to note that, for example, the calculated scatter plots for <7g D 
(Fig. 5) are obtained without any assumption of o"g D = (and cr^ = 0), but the experimental 
exclusion curves for <7g D were traditionally extracted from the data with the spin-neutron (and 
scalar) contribution fully neglected, i.e. under the assumption that <jg D = (and o"gj = 0). This 
one-spin-coupling dominance scheme (always used before the new approaches were pro- 
posed in [115] and in [44, 46, 47]) gave a bit too pessimistic exclusion curves, but allowed direct 
and correct comparison of sensitivities for different DM search experiments. More stringent 
constraints on <7g D can be obtained [44, 46, 47, 115] by assuming both <7g D ^ and <7g D ^ 
(although the contribution of the neutron spin is usually very small because (S£) < (S^)). 
Therefore a direct comparison of an old-fashioned exclusion curve with a new one could in 
principle mislead one to a wrong conclusion about better sensitivity of the new experiment. 

The same conclusion on the one-coupling dominance approach to a great extent concerns 
[46, 47, 51, 52] the direct comparison of the old SI exclusion curves (obtained with a SD = 0) 
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with the new SI exclusion curves (obtained with <t sd > 0) as well as with the results of the 
SUSY calculations. One can see from Fig. 7 that the new-type DAMA/Nal open contour (when 
<t sd > 0) is in agreement with the best exclusion curves of the CDMS and EDELWEISS as 
well as with SUSY calculations. One knows that both of these latter experiments have natural 
germanium (almost purely spinless) as a target and therefore have only little sensitivity to the 
spin-dependent WIMP-nucleon couplings (for them cr gD ~ 0). Therefore, these experiments 
exclude only the pure SI interpretation of the DAMA annual modulation signal [53, 54, 137, 
139, 140]. The statement that this DAMA result is completely excluded by the results of these 
cryogenic experiments and is inconsistent with the SUSY interpretation (see, for example, [141]) 
is simply wrong (see also discussions in [54, 56, 58, 59]). 

The event-by-event CDMS and EDELWEISS background discrimination (via simultaneous 
charge and phonon signal measurements) is certainly very important. Nevertheless the DAMA 
annual signal modulation is one of a few available positive signatures of WIMP-nucleus interac- 
tions and the importance of its observation goes far beyond a "simple" background reduction. 
Therefore, to completely exclude the DAMA result, a new experiment, being sensitive to the 
modulation signal, would have to confirm or exclude this modulation signal on the basis of the 
same or much better statistics. 

Furthermore, taking seriously the positive DAMA result together with the negative results 
of the CDMS and EDELWEISS as well as the results of [57] one can arrive at a conclusion 
about simultaneous existence and importance of both SD and SI WIMP-nucleus interactions. 

V. MIXED SPIN-SCALAR WIMP-NUCLEON INTERACTIONS 

The accurate calculations of spin nuclear structure [31, 66, 70, 88, 98, 100, 101, 106, 108, 112- 
114] demonstrate that contrary to the simplified odd-group approach both (S^ ) and (S^) differ 
from zero, but nevertheless one of these spin quantities always dominates ((S^) <C (S^), or 
(S^) <C (Sp)). It follows form Eq. (6), that if together with the dominance like |(Sp/ n ))| <C 
l(^n(p))l one wou ld have the WIMP-proton and WIMP-neutron couplings of the same order of 
magnitude {not |a„( p )| -C |a p ( n )|), the situation could look like that in the odd-group model and 
one could safely (at the current level of accuracy) neglect a sub-dominant spin contribution in 
the data analysis due to the inequality: \a p (S^)\ <C |a n (S^)|. Nevertheless it was shown in 
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[115] that in the general SUSY model one can meet right a case when a n ( p ) <C a p(n ) and the 
proton and neutron spin contributions are strongly mixed, i.e. \a p (Sp)\ ~ |a n (S^)|. 

To separately constrain the SD proton and neutron contributions at least two new approaches 
appeared in the literature [45, 115]. As the authors of [115] claimed, their method has the 
advantage that the limits on individual WIMP-proton and WIMP-neutron SD cross sections 
for a given WIMP mass can be combined to give a model-independent limit on the properties of 
WIMP scattering from both protons and neutrons in the target nucleus. The method relies on 
the assumption that the WIMP-nuclear SD cross section can be presented in the form <7q n (0) = 
(V^U ± ^) 2 , where M A ar,d ^ are a^Har, q uan titi es, no, directly Jneied 
with measurements. Furthermore, to extract a constraint on the sub- dominant WIMP-proton 
spin contribution one should assume the proton contribution dominance for a nucleus whose spin 
is almost completely determined by the neutrons. From one side, this may look almost useless, 
especially because these sub-dominant constraints are always much weaker than the relevant 
constraints obtained directly with a proton-odd group target (one can compare, for example, 
the restrictive potential of the NAIAD exclusion curves in Figs. 5 and 6). From another side, the 
very large and very small ratios u p /a n ~ |a p |/|a n | obtained in [115] correspond to neutralinos 
which are extremely pure gauginos. In this case Z-boson exchange in SD interactions is absent 
and only sfermions give contributions to the SD cross sections. This is a very particular (fine- 
tuning) case which is hardly to be in agreement with the present SUSY search experiments. 
Following an analogy between neutrinos and neutralinos one could assume that neutralino 
couplings with the neutron and the proton should not be very different and one could expect 
preferably |a n |/|a p | ~ 0(1). The relation |a n |/|a p | pa 0(1) was checked in [69, 121] for large 
LSP masses. For relatively low LSP masses m x < 200 GeV in effMSSM [8, 37, 116-119, 124] 
the a n -to-a p ratio is located within the bounds [51]: 

a, 



0.5 < 



a p 



< 0.8. (23) 



Therefore in the model the couplings are almost the same and one can safely neglect the (Sp 4 ^)- 
spin contribution in the analysis of the DM data for a nuclear target with (S^) <C {^n( P ))- 

Furthermore, when one compares in the same figure an exclusion curve for SD WIMP- 
proton coupling obtained without sub-dominant SD WIMP-neutron contribution and without 
SI contribution (all curves in Fig. 5 except the one for NAIAD [136] and the one for Tokyo-LiF 
[142]), with a curve from the approach of [115], when the sub-dominant contribution is included 
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(the NAIAD and Tokyo-LiF curves in Fig. 5), one "artificially" improves the sensitivity of the 
latter curves (NAIAD or Tokyo-LiF) in comparison with the former ones. To be consistent and 
for reliable comparison of sensitivities of these experiments, one should, at least, coherently 
recalculate all previous curves in the new manner. This message was clearly stressed in [46]. 

The same arguments are true for the results of the SIMPLE experiment [143] and search 
for DM with NaF bolometers [144] where the SI contribution seems also to be completely 
ignored. Although 19 F has the best properties for investigation of WIMP-nucleon spin- 
dependent interactions (see, for example [88]) it is not obvious that one should completely 
ignore spin-independent WIMP coupling with the fluorine. For example, in the relation 

which follows from Eqs. (4)-(6), it is not a priori clear that 



~ u 



SD 

A,n 



4 + + 



— f^F <C i-e. the SI WIMP-nucleus interaction is much weaker than the sub-dominant SD 

°SD CT SD 

WIMP-nucleus one. At least for isotopes with atomic number A > 50 [29, 33] the neglection 
of the SI contribution would be a larger mistake than the neglection of the sub-dominant SD 
WIMP-neutron contribution, when the SD WIMP-proton interaction dominates. 

Therefore we would like to note that the "old" odd-group-based approach to analysis of 
the SD data from experiments with heavy enough targets (for example, Ge-73) is still quite 
suitable, especially when it is not obvious that (both) spin couplings dominate over the scalar 
one. 

From measurements with Ge one can extract, in principle, not only the dominant constraint 
for WIMP-nucleon coupling a n (or <7g D ) but also the constraint for the sub-dominant WIMP- 
proton coupling a, p (or Og D ) using the approach of [115]. Nevertheless, the latter constraint will 
be much weaker in comparison with the constraints from p-odd group nuclear targets, like 19 F 
or I. This fact is illustrated by the NAIAD (Nal, 2003) curve in Fig. 6, which corresponds to 
the sub-dominant WIMP-neutron spin contribution extracted from the p-odd nucleus 127 1. 

Another approach for the mixed spin-scalar coupling data presentation, of Bernabei et al. 
[45] , is based on an introduction of the so-called effective SD nucleon cross section Og^ (0) (ctsd 
in [45, 46]) and coupling mixing angle 9 (see Eq. (7)) instead of 0^(0) and ag D (0). With these 
definitions the SD WIMP-proton and WIMP-neutron cross sections are given by relations (8). 

In Fig. 8 the WIMP-nucleon spin and scalar mixed couplings allowed by the annual modu- 
lation signature from the 100-kg DAMA/Nal experiment are shown inside the shaded regions. 
The regions from [45, 46] in the (<7si, 0sd) space for 40 GeV< m WIMP <110 GeV cover spin- 
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scalar mixing coupling for the proton (9 = case of [45, 46], left panel) and spin-scalar mixing 
coupling for the neutron (9 = 7r/2, right panel). From nuclear physics one has for the proton 
spin dominated 23 Na and 127 I ||4 < 0.1 and |^ < 0.02 + 0.23, respectively. For 9 = due to 
the p-oddness of the I target, the DAMA WIMP-proton spin constraint is the most severe one 
(see Fig. 5) . In the right panel of Fig. 8 we also present the exclusion curve (dashed line) for 
the WIMP-proton spin coupling from the proton-odd isotope 129 Xe obtained under the mixed 
coupling assumptions [45] from the DAMA-LiXe (1998) experiment [145-147]. For the DAMA 
Nal detector the 9 = n/2 means no (S p ) contribution at all. Therefore, in this case DAMA gives 
the sub-dominant (S n ) contribution only, which could be compared further with the dominant 
(S„) contribution in 73 Ge. The scatter plots in Fig. 8 give Ugj as a function of crg D (left panel) 
and o"g D (right panel) calculated in the effMSSM with parameters from Eq. (22) under the same 
constraints on the relic neutralino density as in Figs. 4 [51]. Filled circles (green) correspond 
to relic neutralino density 0.0 < Q x ti^ < 1.0, squares (red) correspond to sub-dominant relic 
neutralino contribution 0.002 < £\/iq < 0.1 and triangles (black) correspond to the WMAP 
density constraint 0.094 < Q x h 2 < 0.129 [18, 19]. 
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FIG. 8: The DAMA-Nal allowed region from the WIMP annual modulation signature in the (£osi, 
£csd) space for 40 < rn wlMP < 110 GeV [45, 46]. The left panel corresponds to the dominating (in 
127 I) SD WIMP-proton coupling alone [6 = 0) and the right panel corresponds to the subdominating 
SD WIMP-neutron coupling alone (6 = it/ 2). The scatter plots give correlations between Ogj and cr SD 
in the effMSSM (£ = 1 is assumed) for m x < 200 GeV [51]. In the right panel also the DAMA liquid 
xenon exclusion curve from [45] is given (dashed line). From [51]. 
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The constraints on the SUSY parameter space in the mixed coupling framework in Fig. 8 
look, in general, much stronger in comparison with the traditional approach based on the 
one-coupling dominance. It follows from Fig. 8 that when the LSP is the subdominant DM 
particle (squares in the figure), SD WIMP-proton and WIMP-neutron cross sections at a level 
of 3 -T- 5 • 1CT 3 pb are allowed, but the WMAP relic density constraint (triangles) together with 
the DAMA restrictions leaves only < 3 • 1CT 5 pb without any visible reduction of allowed 
values for Ogj. In general, according to the DAMA restrictions, very small SI cross sections are 
completely excluded, only <7gj > 3 4- 5 ■ 1CT 7 pb are allowed. The SD cross section is not yet 
restricted at all. It is seen that for the allowed values of the SI contribution the SD DAMA 
sensitivity did not yet reach the upper bound for the SD LSP-proton cross section of 5 • 10~ 2 pb 
calculated for the nucleon spin structure from [148]. 

In general, the famous DAMA "conflict" with the other (negative) DM results can be safely 
bypassed on the basis of the above-mentioned mixed spin-scalar coupling approach, where both 
SD and SI couplings are considered simultaneously as non- negligible. 

VI. THE MIXED COUPLING APPROACH FOR THE HIGH-SPIN 73 Ge 

In this section, on the basis of the high-spin 73 Ge detector HDMS [91, 92], the mixed spin- 
scalar coupling approach is used to demonstrate how one can significantly improve the quality 
of the exclusion curves in comparison with the one-coupling dominance result of [79]. 

The Heidelberg Dark Matter Search (HDMS) experiment used a special configuration of two 
Ge detectors to efficiently reduce the background (due to anti-coincidence of inner and outer 
detectors) [79, 91, 92]. A small, p-type Ge crystal (enriched by 86% in 73 Ge) is surrounded by 
a well-type natural Ge crystal, both being mounted into a common cryostat system (see left 
panel in Fig. 9 for a schematic view). The HDMS with enriched 73 Ge inner detector was the 
first and till now unique setup with a high-spin (J=9/2) Ge target isotope for direct DM search. 
The main idea of the new combined analysis relies on the unique possibility that two different 
isotope targets (from natural Ge and enriched 73 Ge) were used as inner detector in the same 
HDMS setup under the same outer background conditions of LNGS. 

In fact the first simple estimation of the prospects for DM search and SUSY constraints with 
the high-spin 73 Ge detector HDMS assuming mixing of WIMP-neutron spin and WIMP-nucleon 
scalar couplings together with available results from the DAMA-Nal and LiXe experiments 
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FIG. 9: Left: Schematic view of the HDMS detector configuration. A small Ge crystal is surrounded 
by a well type Ge-crystal, the anti-coincidence between them is used to suppress background created 
by external photons. Right: The HDMS detector in its open shield during the installation in the Gran 
Sasso Underground Laboratory. The inner shield is made of 10 cm of electrolytic copper, the outer 
one of 20 cm of Boliden lead. From [92, 149]. 



[44, 46, 47, 145-147] was performed in [51]. Furthermore, recently in the mixed spin-scalar 
coupling approach the data from both HDMS experiments with natural Ge and with enriched 
73 Ge were simultaneously re-analyzed. This new analysis together with a new procedure for 
background identification and subtraction from the measured 73 Ge spectrum allowed one to 
obtain a significant (about one order of magnitude) improvement for the limits on the WIMP- 
neutron spin-dependent coupling. As a result the HDMS experiment is now giving the most 
sensitive limits on the WIMP-neutron spin coupling for WIMP masses larger than 60-65 GeV/c 2 
[93]. 

The evaluation of the DM limits (exclusion curves) on the WIMP-nucleon SD or SI cross 
section follows, in general (see, for example [91, 150, 151]), the conservative assumption that 
the whole experimentally measured spectrum is saturated by the WIMP-induced events. Con- 
sequently, any excess events from the calculated spectrum above the relevant experimental 
spectra in any energy interval are considered as forbidden (at a given confidence level). In our 
case we assume that for any given WIMP mass m x both crgi = o~si( m x) an d °"sd = ctsd(^x) 
WIMP-nucleon interaction cross sections are excluded if 



dR( E ) dR 1 ^ 3 ^ 8 " 
"' ' U = K sl (E R ,m x )asi + K SB (E R ,m x )a S B > (e, Er) = # data (e, E ¥ 
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or both upper limit values for <7si and <tsd can be obtained as solutions of the following equation 



k si (E r , m x ) a S i + k sd (E r , m x ) a SD = i? data (e, E R 



for all available recoil energies E R above the threshold energy e. The notations are given by 
Eqs. (8), (7), (15), and (16). The sub-dominant contribution from WIMP-proton spin coupling 
proportional to (S^ ) can be safely neglected for 73 Ge. The 73 Ge isotope looks with a good 
accuracy like an almost pure neutron-odd group nucleus with (S„) (S p ) (Table II). Therefore 
in our consideration a SD = o"g D and cos# = 0. For the WIMP mass density in our Galaxy the 
value p x = 0.3 GeV/cm 3 is used. 

To find both a sl and <t sd , for any given m x , in accordance with Eq. (24) the following 
functional 

2 

(25) 
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can be numerically minimized, where R?(Ei) and ARi(Ei) are measured rate and its error (in 
counts/day/kg/keV) in z-th energy bin for j-th used spectrum (j = 1,2 for natural Ge and 
Ge-73). The two main used spectra are given in Fig. 10. Only the "cleanest" background 
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FIG. 10: HDMS spectra used in the analysis. Left panel: Spectrum obtained in the first prototype 
phase of the HDMS experiment with the inner detector from natural Ge (exposure 9.9 kg d) [150]. 
Right panel: Spectrum obtained from the full HDMS setup with the inner detector from enriched 73 Ge 
(exposure 85.5 kg d) [79]. It is separated in 3 sets of runs of 30.9, 29.5 and 27.6 kg d, respectively. 
The blue histogram corresponds to the latest (most clean) runs 721-1000. From [79]. 



spectrum with the Ge target collected in the latest runs of the experiment (with numbers 
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721-1000) was used in the analysis (blue color histogram in the right panel of Fig. 10). For 
both spectra the visual energy threshold e = 4 keV is used. 

Comparing both the 'most accurate' (blue, runs 721-1000 in Fig. 10) HDMS 73 Ge-spectrum 
and the natural Ge spectrum one can see obviously some non-vanishing extra background 
contribution in the first spectrum relative to the second one. In general, such a possibility is 
not new. The improvement in the exclusion curves by taking into account known sources of 
background during DM searches with Ge detectors was demonstrated, for example, in [152] and 
further discussed in [153]. Therefore for this HDMS 73 Ge-spectrum we allow the possibility to 
fit simultaneously with the SD and SI cross sections some constant (as function of the recoil 
energy) background contribution, too. The effect of this extra background contribution is 
discussed later on. 

For a semiconductor germanium detector one has to take into account the ionization effi- 
ciency For the HDMS Ge setup in the (visual) energy interval 4 < E R < 50 keV a simple rela- 
tion between the visible recoil energy and the real recoil energy — E vis = E R = 0.14_E r 1 c 3 tl 9 rccoil 
QE veal recoil — can be used with Q = 0.33 being the quenching factor for Ge [79, 91, 150, 151]. 

One can note that for any WIMP mass, m x , and any target mass, Ma, due to kinematics 
one has not to expect any WIMP-induced event at all with 

rp ^ p max _ n 2t, max^ m x _ n 2 "L 2 „ _ „ 

> -&R — 72 — V^TT — V-Ai ^max ~ ^esc- I/O J 

[Ma + m x ) Ma 

For example, £ R max = 4 (50) keV, for a Ge detector and m x = 12 (77) GeV/c 2 . It is clear from 
Eq. (26) that for fixed v max , Ma and the detector energy threshold e there are undetectable 
WIMPs (with rather light m x ) if the maximal recoil energy, they can produce, is smaller than 
the threshold: E R max < e. Therefore one has two restrictions for the theoretical event rate as 
function of the WIMP mass m x : 

R(E R > £ R max (m x )) = and R(E R max (m x ) < e) = 0. (27) 

The first one allows background (for the WIMP-nucleus signal) estimation in the 73 Ge measured 
spectrum, which could lead to a remarkable improvement of the deduced exclusion curves. 

Now we turn to our main analysis of both HDMS spectra in the mixed spin-scalar coupling 
approach and extract limits for both cross sections <r SI and <r SD simultaneously using formulas 
(24)-(25). To obtain from the available data the most accurate exclusion curve one can use 
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two minimization approaches. The first (main) approach relied on indeed direct simultaneous 
determination of the SD and SI WIMP-nucleon limits for a given WIMP mass m x (exclusion 
curves for c SD (m x ) and c si (m x )) by means of minimization of the discrepancy between our 
calculated estimations of the expected rates and both above-mentioned experimental spectra. 
In our second (auxiliary) approach assuming SI coupling dominance (<r SD = 0) we first extracted 
only the SI WIMP-nucleon limit cr gI (m x ) from the natural Ge spectrum (left panel in Fig. 10). 
Next, for each WIMP mass and above defined <J SI (m x ) we extracted only the SD WIMP-nucleon 
limit o" S i( m x) fr° m the cleanest background spectrum with 73 Ge (Runs 721-1000, 27.6 kg d, 
blue spectrum in right panel of Fig. 10). In both cases the obtained results are rather similar. 

Furthermore trying to improve the quality of the exclusion curves one can use a sliding 
variable energy window to check the excess events above the experimental spectrum (in these 
energy window intervals) as used in previous papers [79, 91, 92, 150]. The minimum among 
the cross section values obtained via the multiple fits is taken as the cross section for the 
corresponding WIMP mass. We used 5 keV minimal width of this energy window as in [150, 151] 
and a 10-keV window as well. 

First, the possible improvements of the exclusion curves due to variation of minimization 
conditions, used in the data analysis, were studied. The relevant exclusion curves obtained 
from the simultaneous analysis of both HDMS spectra within the mixed spin-scalar coupling 
approach are given in Fig. 11 as function of the WIMP mass. 

The left panel of Fig. 11 shows the upper limits for the spin-independent WIMP-nucleon 
cross section cxsi as function of WIMP mass obtained under different minimization conditions. 
The right panel of Fig. 11 shows the upper limits for the spin-dependent WIMP-neutron cross 
section <7sd (in our approximation <t sd = o"sd) as function of WIMP mass which correspond 
to the above-mentioned <tsi limits from the left panel. For example, the (red) exclusion curve 
labeled with "si-sd-05(l-5)" presents in the left panel the SI limits o" S i( m x) extracted only 
from the natural Ge spectrum (left panel, Fig. 10) with 5-keV sliding windows within only the 
first 5 lowest energy windows and under the assumption <r SD = 0. The relevant (red) "si-sd- 
05(1-5)" exclusion curve in the right panel shows correlated SD limits cr SD (m x ) extracted from 
the best enriched 73 Ge spectrum (Runs 721-1000, blue spectrum, right panel, Fig. 10) when 
cr SI = o"si( m x) 7^ fro m the left panel. The (thick blue) curves labeled with "sd+si-05(l-5)" are 
obtained from simultaneous minimization of both above-mentioned spectra from Fig. 10 using 
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FIG. 11: Spin- independent <tsi (left panel) and spin-dependent <t sd (right panel) cross section upper 
limits as function of WIMP mass (exclusion curves) obtained from simultaneous analysis of both 
HDMS Ge spectra under different minimization conditions. The label "si-sd-05(l-5)" means that 
only the spectrum from natural Ge (left panel, Fig. 10) is used first to extract a SI limit (auxiliary 
approach) using 5-keV sliding window and taking into account only the first 5 lowest energy windows 
to obtain the limits. Curves (blue) labeled with "sd+si-05(l-5)" are obtained from indeed simultaneous 
minimization of both natural Ge and the best (right panel, Fig. 10) enriched Ge spectra using the 
5-keV sliding window and only the first 5 lowest energy windows. The label "si+sd-10(2-3)" denotes 
the same procedure, but with 10-keV sliding window and with the 2nd and the 3rd lowest energy 
windows. The other labels have analogical meaning. Dashed lines correspond to some other analyses 
and are given for comparison ("HDMS2000, nat Ge" from [150], "HDMS05/SD" from [79]). The 
thin exclusion curves "sd-45", "sd-10(l-8)" and "sd-05(l-10)" are obtained from the traditional one- 
coupling dominance fit of the enriched Ge spectrum only and are given to illustrate the role of the 
sliding window width as well as the consistency with the previous result of [79]. 



the 5-keV sliding window and the first 5 lowest energy windows. The other exclusion curves 
(with labels "sd+si-10(2-3)" and "sd-si-05(l-3)") are given in Fig. 11 to illustrate the exclusion 
curve dependence on the width of the sliding window and optimal (or non-optimal) choice of 
minimization regions. 

The dashed lines correspond to some other analyses and are given for comparison 
("HDMS2000, nat Ge" from [150], "HDMS05/SD" from [79]). The thin exclusion curves "sd- 
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45", "sd-10(l-8)" and "sd-05(l-10)" are obtained from the traditional one-coupling dominance 
fit of the enriched 73 Ge spectrum only and are given also to illustrate the role of the width 
of the sliding window as well as the consistency with the previous result of [79]. The small 
discrepancy between the "HDMS05/SD" and "sd-05(l-10)" curves for low WIMP masses is 
mainly due to our restrictions (27). 

The pair of SI and SD exclusion curves with label "sd+si-05(l-5)" corresponds to the optimal 
parameters of the fitting procedure and presents the best correlated exclusion curves obtained 
from simultaneous minimization of both spectra from natural Ge spectrum and the "cleanest" 
enriched 73 Ge spectrum. The pair has the best SI exclusion curve (left panel) simultaneously 
with almost the best SD exclusion curve (right panel). The visible difference between the best 
SI curve "sd+si-05(l-5)" and the "HDMS2000, nat Ge" curve from [150] is due to restrictions 
(27) and another energy threshold (see below) used in our analysis. 

Therefore, from the right panel of Fig. 11 one can conclude that the most sensitive exclusion 
curves for the WIMP-neutron spin interaction ("si-sd-05(l-3)", black and "si+sd-05(l-5)", blue) 
improve the relevant one-coupling dominance limit of "HDMS05/SD" [79] within a factor of 
2-3 depending on the WIMP mass. This is a clear result of the mixed spin-scalar approach. 

Now we consider two other possibilities to improve the quality of the exclusion curves ex- 
tracted from both Ge spectra within the mixed coupling scheme. The first one is a lower recoil 
(visible) energy threshold for the natural Ge detector of HDMS setup. The second one is a new 
procedure of background subtraction from the measured spectrum of the 73 Ge isotope. This 
procedure strongly relies on the existence of a really measured spectrum. 

In Figure 12 the thin (black) exclusion curve labeled with "SD only fit" is the result of the 
one coupling dominance analysis of the enriched 73 Ge spectrum only. It repeats the relevant 
curve (labeled with "sd-05(l-10)") from Fig. 11. As mentioned before, the curve is consistent 
with the previous analysis of [79], given here with curve "HDMS05/SD". 

The curves labeled with "no bg subtraction" repeat here the best exclusion curves ( "si+sd- 
05(1-5)" in Fig. 11) obtained within the mixed spin-scalar coupling approach from simultaneous 
analysis of both natural Ge and the best enriched 73 Ge spectra. The same recoil energy thresh- 
old of 4 keV was taken for both Ge spectra. This threshold corresponds to the real threshold 
of the HDMS final setup with enriched 73 Ge [91]. 

We reproduce in Fig. 12 both exclusion curves "SD only fit" and "no bg subtraction" (from 
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FIG. 12: Exclusion curves from simultaneous fit of the data from both HDMS setups. Left panel: 
Spin-independent cross section upper limits <tsi in pb as function of WIMP mass in GeV. Right panel: 
Spin-dependent cross section upper limits ctsd as function of WIMP mass which correspond to the 
ctsi limits from the left panel. The label "no bg subtraction" shows the best exclusion curves (in red) 
obtained from simultaneous minimization of both natural Ge (left panel, Fig. 10) and the best (Runs 
721-1000, blue spectrum, right panel, Fig. 10) enriched 73 Ge spectra within mixed spin-scalar coupling 
approach without any background subtraction from the Ge spectra. Curves (in blue) labeled with "4 
keV threshold" are from the same simultaneous minimization of both natural Ge and the best enriched 
73 Ge spectra but with extra background of 0.11 events/kg/day extracted from the 73 Ge spectrum. The 
curves (in black) labeled with "2 keV threshold for nat Ge" are from the same procedure as above, 
but when threshold for the natural Ge spectrum is equal to 2 keV. The curve reproduces the best 
HDMS SI limits from [150] given here as "HDMS2000, nat Ge" (dashed green). The thin (black) 
exclusion curve "SD only fit" from the traditional one-coupling dominance fit of the enriched 73 Ge 
spectrum only is given for comparison with previous analysis of [79] labeled with "HDMS05/SD". The 
best exclusion curve for WIMP-neutron spin coupling from CDMS collaboration [154] is labeled with 
"CDMS-2006/SD" (dash green). The last dashed (green) line "DAMA98/Xe" corresponds to DAMA 
results from [146]. Another comparative result from ZEPLIN-I [155] (not shown) is located above the 
CDMS curve nearby the black one. 

Fig. 11) for our further consideration and with the aim to clearly demonstrate again that 
the most sensitive HDMS exclusion curve ("no bg subtraction") for the WIMP-neutron spin 
interaction improves the relevant one-coupling dominance limit of "HDMS05/SD" [79] and the 
"SD only fit" curve within a factor of 2-3. 
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We stress that this "no bg subtraction" curve is obtained from the raw HDMS data without 
any active or passive background substraction. 

The visible difference at low WIMP masses between SI exclusion curve "no bg subtraction" 
and the "HDMS2000, nat Ge" curve from [150] (in the left panel of Fig. 12) is mainly due to 
the lower recoil energy threshold of 2 keV used in [150] for the natural Ge detector. The curve 
"2 keV threshold for nat Ge" (black in the left panel of Fig. 12) obtained indeed with a 2 keV 
energy threshold for the spectrum of natural Ge proves the reason of the difference. 

Furthermore, the real measured spectrum of enriched 73 Ge (Runs 721-1000, blue spectrum, 
right panel, Fig. 10) and the first relation (27) allow one to estimate some number of counts 
in the spectrum which can not be produced by means of any WIMP-nucleus interaction. In 
accordance with (27) for any m x , fixed f max and Ma there is a maximal recoil energy E'p^ ax {rn-^) 
(26) for which WIMP-nucleus interactions are unable to produce any signal if En > i?R max (m x ) 
(i.e. when the measured recoil energy is larger than the maximally possible recoil energy for 
a given WIMP mass). Therefore, for fixed m x the measured recoil spectrum in the region 
E K > E R max (m x ) is directly some background which can be approximated, for example, as a 
constant function of the recoil energy, independent of m x . One can estimate these background 
constants for each allowed m x (still £ , R max (m x ) < 50 keV) and assume the minimal of these 
constants (0.11 events/kg/day/keV) to be the mean background for all measured Er and all 
m x . Physically this extra background is completely independent on any WIMPs, therefore 
being estimated for rather small m x < 100 GeV, it can be used for all WIMP masses as well. 

Therefore, with common energy threshold of 4 keV, the simultaneous minimization of both 
natural Ge and the best 73 Ge spectra with the above-mentioned extra background of 0.11 
events/kg/day has supplied us with the pair of SD and SI exclusion curves, labeled with "4 
keV threshold" in Fig. 12. This SD curve improves (at least for m x > 60 GeV) currently the 
best exclusion curve (labeled with "CDMS-2006/SD" , dash green) for the WIMP-neutron spin 
coupling from the CDMS collaboration [154]. The other dashed lines correspond to some other 
analyses and are given for comparison ("HDMS2000, nat Ge" from [150], "HDMS05/SD" from 
[79], and "DAMA98/Xe" from [146]). 

The curves (in black, right panel of Fig. 12) labeled with "2 keV threshold for nat Ge" are 
from the same fit procedure with the extra background, but when the threshold for the natural 
Ge spectrum is equal to 2 keV. In this mentioned above, one reproduces the best 
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HDMS SI limits "HDMS2000, nat Ge" (dashed green) from [150]. As it is seen, the background 
subtraction from the 73 Ge spectrum only very weakly affects the correspondent SI curves in 
Fig. 12 (left panel). 

The main results of the analysis performed in the mixed spin-scalar coupling approach are 
the (correlated) limits for the cross sections a SI and cr SD . Indeed, despite the traditional form 
of presentation of the SD and SI exclusion curves in Fig. 12 as function of WIMP mass one 
should keep in mind that these a sl and <7g D constraints for fixed WIMP mass are strongly 
correlated. This correlation is presented explicitly in Fig. 13, where the dependence on WIMP 
mass is given indirectly by means of the points running over the curves and marked with 
relevant WIMP mass values. For example, the point marked with label "200 GeV" shows 
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FIG. 13: Correlated spin-dependent cross section upper limit ctsd and spin-independent cross section 
upper limit <tsi obtained from the simultaneous analysis of the data from both HDMS setups. Points 
running over the curves mark the relevant WIMP mass values. For example, the point marked with 
label "200 GeV" shows the best simultaneous upper limits for SD and SI WIMP-nucleon interaction 
cross sections for WIMP mass 200 GeV/c 2 . For the first time similar plots were given in [45]. From 
[93]. 
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the simultaneous upper limits for SI and SD WIMP-nucleon interaction cross sections <r SI and 
<7g D for m x = 200 GeV/c 2 . This, in principle, gives one a new requirement (for a SUSY- 
like theory) that for any fixed WIMP mass m x one should have <r SI (theor.) < a sl (fitted) and 
(Tg D (theor.) < (Tg D (fitted) simultaneously. 

For the sake of completeness one can compare the limits obtained (in the spin-scalar mixed 
coupling approach) from the HDMS experiments on the SD and SI WIMP-nucleon interaction 
with the relevant constraints extracted by the DAMA collaboration from measurement of the 
annual signal modulation with Nal target [46]. Following the DAMA positive evidence one can 
accept that the most prefered interval of the WIMP mass is 40 GeV/c 2 < Tn wlMP < HO GeV/c 2 . 
Other consequences of the fact one can find in [51, 52]. 

In Fig. 14 the new HDMS-2006 limits on the SD and SI WIMP-nucleon interactions are 
compared with the relevant DAMA constraints extracted from measurement of the annual 
signal modulation with a Nal target [46] as well as with calculation in the low-energy effective 
MSSM [51]. To perform the comparison with the DAMA allowed region of a sl and <Tg D (Fig. 14) 
only one WIMP mass m x = 80 GeV/c 2 (the star in the figure) was chosen for illustration. The 
upper limits for all other WIMP masses, 40 GeV/c 2 < m x < 110 GeV/c 2 , are very close to 
this point. The point marked with the star gives our simultaneous upper limits for (T SI and 
<7g D for m x = 80 GeV/c 2 . For this WIMP mass, values of a sl located above the horizontal line 
and of <Tg D located to the right side of the vertical line are excluded. Therefore, for any fixed 
WIMP mass in the domain 40 GeV/c 2 < m x < 110 GeV/c 2 SUSY-like calculations should 
give simultaneously a sl below the relevant (to the fixed WIMP mass) horizontal line and <7g D 
to the left side of the relevant vertical line, respectively. These limits improve the DAMA-Xe 
limit significantly (about one order of magnitude) and exclude a DAMA allowed region of large 
spin-dependent WIMP-neutron cross sections. 

Here, perhaps, is the right place to make another remark concerning the possibility to com- 
pare results from DM search experiments with passive background reduction (like DAMA, 
HDMS, etc) and experiments with (mostly) active background reduction (like CDMS, EDEL- 
WEISS, etc). First, we note, that obviously any extra positively defined background-like con- 
tribution to the spectra will decrease the extracted (upper limit) values of the SD and SI cross 
section. Next, within the passive background reduction scheme the measured spectrum is not 
affected by hard- or software influence during the data taking. Extra further background re- 
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FIG. 14: The DAMA-Nal allowed region (inside yellow band) for SD WIMP-neutron coupling versus 
SI WIMP-nucleon coupling is from [46] and corresponds to 40 GeV/c 2 < m x < 110 GeV/c 2 . The 
scatter plots from [51] give correlations between Ogj and (Xg D in the effMSSM for m x < 200 GeV. 
The squares (red) correspond to sub-dominant relic neutralino contribution 0.002 < f2y/ig < 0.1 and 
triangles (black) correspond to WMAP relic neutralino density 0.094 < O x /ig < 0.129. The dashed 
line from [45] shows the DAMA-LiXe (1998) exclusion curve for m x = 50 GeV/c 2 . The star gives 
our simultaneous upper limits from the HDMS experiment [93] for <t si and <7g D for m x = 80 GeV/c 2 . 
Therefore values of <t si above the horizontal line and of <7g D located right from the vertical line are 
excluded by our analysis. 

duction can be done off-line on the basis of careful investigation of the spectrum itself or, for 
example, with help of pulse shape analysis. In this case the extracted background contribution 
is under control and well denned. On the other side, within the active background reduction 
approach the measured spectrum already contains results of this active reduction influence on 
the data taking process. In this case it is not simple to hold under control the real level of 
extracted on-line background contribution which easily can be overestimated (see, for exam- 
ple, the recent discussion of the ZEPLIN-I sensitivity in [156]). Therefore, due to this obvious 
difference a direct comparison of exclusion curves from experiments with passive and active 
background reductions could be, in principle, rather misleading. 
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VII. DISCUSSIONS 

The problem of the dark matter in the Universe is a challenge for modern physics and ex- 
perimental technology. To solve the problem, i.e. at least to detect dark matter particles, 
one simultaneously needs to apply the front-end knowledge of modern Particle Physics, Astro- 
physics, Cosmology and Nuclear Physics as well as one should develop, and use over long time 
extremely high-sensitive experimental setups, and complex data analysis methods. 

Weakly interacting massive particles (WIMPs) nowadays are among the best motivated non- 
baryonic dark matter candidates. In particular, the lightest neutral supersymmetric particle 
(LSP), the neutralino, is a very good WIMP candidate. The motivation for supersymmetry 
arises naturally in modern theories of particle physics. 

To estimate the expected direct detection rate for these WIMPs any SUSY-like model, for 
example, an effective low-energy minimal supersymmetric extension of the Standard Model 
(effMSSM), or some measured data, for example, from the DAMA experiment [157], can be 
used. On this basis the WIMP-proton and WIMP-neutron spin and scalar cross sections at 
zero-momentum transfer (og^ (0) and <7gj n (0)) can be calculated. These calculations one usually 
compares with measurements, which (with the only exception of the DAMA result) are pre- 
sented in the form of exclusion curves — upper limits of cross section as functions of the WIMP 
mass. In the case of non-observation of any DM signal an exclusion curve simply reflects the 
sensitivity of a given direct DM search experiment and potentially allows one to constrain some 
version of SUSY-like theory, if the curve is sensitive enough. Therefore the best exclusion curve 
is currently a clear aim of almost all dark matter search experiments (DAMA, LIBRA, and GE- 
NIUS perhaps are/were the only exceptions). The main competition between the experiments 
runs in the field of these exclusion curves. 

Before 2000 all exclusion curves were evaluated mainly in the one-coupling dominance ap- 
proach (when only one cross section was defined from the measured spectra for fixed WIMP 
mass), which gave slightly pessimistic (for spin- non- zero target experiments), but universal 
limits for all experiments. One would say that the competition between DM experiments was 
honest. The predictions from SUSY-like models were in general far from being reached by the 
data. 

Mainly after the paper [115] was published in 2000 (and as well after the DAMA evidence 
[46]) a new kind of exclusion curves appeared. In particular, for the first time these curves 
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were obtained for the spin-dependent WIMP-nucleon cross section limits when non-zero sub- 
dominant spin WIMP-nucleon contributions were taken into account [136, 142]. This procedure 
obviously improved the quality of the exclusion curves. Therefore a direct comparison of an 
old-fashioned exclusion curve with a new one could in principle mislead one to a wrong conclu- 
sion about better sensitivity of the more recent experiments. There is generally some possible 
incorrectness in the direct comparison of the exclusion curves for the WIMP-proton(neutron) 
spin-dependent cross section obtained with and without the non-zero WIMP-neutron (proton) 
spin-dependent contribution. Furthermore the above-mentioned incorrectness concerns to a 
great extent the direct comparison of spin-dependent exclusion curves obtained with and with- 
out non-zero spin-independent contributions [46, 47]. Taking into account both spin couplings 
a p and a n but ignoring the scalar coupling cq, one can easily arrive at a misleading conclusion 
especially for not very light target nuclei when it is not obvious that (both) spin couplings 
dominate over the scalar one. To be consistent, one has to use a mixed spin-scalar coupling 
approach as for the first time proposed by the DAMA collaboration [44, 46, 47]. 

It was argued in [44, 46, 47, 51, 52] that potentially misleading discrepancies between the 
results of different dark matter search experiments (for example, DAMA vs CDMS and EDEL- 
WEISS) as well as between the data and the SUSY calculations can be avoided by using the 
mixed spin-scalar coupling approach, where the spin-independent and spin-dependent WIMP- 
nucleon couplings are a priori considered to be both non-zero. 

The mixed spin-scalar coupling approach was applied to analyze the data from both HDMS 
experiments with natural Ge and with the neutron-odd group high-spin isotope 73 Ge. The 
approach allows both upper limits for spin-dependent (Jg^ and spin-independent a sl cross 
sections of WIMP-nucleon interaction to be simultaneously determined from the experimental 
data. In this way visible improvement in form of exclusion curves is achieved relative to the 
traditional one-coupling dominance scheme [93]. The agreement of the obtained <7g D and o~ SI 
with parameter regions allowed from the observation of the annual modulation signature by 
the DAMA collaboration is demonstrated. The above-mentioned correlations between <jg D and 
(J SI can be considered as a new requirement, which demands that for any fixed WIMP mass 
m x one should have a SI (theor.) < <t si (fitted) and <jg D (theor.) < <7g D (fitted), simultaneously, 
provided ag-^gj-^theor.) are calculated in any underlying SUSY-like theory. For the first time 
a similar result for Nal was mentioned by the DAMA collaboration [45]. 

It is important to note, that without proper knowledge of the nuclear and nucleon structure 
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it is not possible to extract reliable and useful information (at least in form of these <7g D 
and er SI cross sections) from direct dark matter search experiments. However, astrophysical 
uncertainties, in particular the DM distribution in vicinity of the Earth [53-59], make the 
problem of interpretation of results of the DM search experiments far more complicated. At the 
moment to have a chance to compare sensitivities of different experiments people adopted one 
common truncated Maxwellian DM particle distribution, but nobody can prove its correctness. 
Only in the case of indeed direct DM detection one can make some conclusions about the real 
DM particle distribution in the vicinity of the Earth. 

Furthermore, almost by definition (from the very beginning), a modern experiment aiming 
at the best exclusion curve is doomed to non-observation of the DM signal. This is due to the 
fact, that a typical expected DM-signal spectrum exponentially drops with recoil energy and it 
is practically impossible to single it out from a background non-WIMP spectrum of a typical 
(semiconductor) detector, which is as usual exponential as well. 

In fact, one needs some clear, or "positive" signature of WIMP particles interactions with 
target nuclei. Only exclusion curves are not enough. Ideally this signature should be a unique 
feature of such an interaction (see for example [158]). 

There are some typical characteristics of WIMP particle interactions with a nuclear target 
which can potentially play the role of such positive WIMP signatures (see for example [159]). 
First of all WIMPs produce nuclear recoils, while most radioactive backgrounds produce elec- 
tron recoils. Nevertheless, for example, neutrons (and any other heavy neutral particle) also 
can produce nuclear recoils. There exist also some proposals which rely on WIMP detection 
via electron recoils (see for example [160-162]). 

Due to the extremely rare event rate of the WIMP-nuclear interactions (the mean free path 
of a WIMP in matter is of the order of a light-year) one can expect two features. The first one 
is that the probability of two consecutive interactions in a single detector or two closely located 
detectors is completely negligible. Multiple interactions of photons, gamma-rays or neutrons 
under the same conditions are much more common. Therefore only non-multiple interaction 
events can pretend to be from WIMPs. The second one is a uniform distribution of the WIMP 
induced events throughout a detector. This feature can also be used in future to identify 
background events (from photons, neutrons, beta and alpha particles) in rather large-volume 
position-sensitive detectors. 
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The shape of the WIMP-induced recoil energy spectrum can be predicted rather accurately 
(for given WIMP mass, for fixed nuclear structure functions and astrophysical parameters). 
The observed energy spectrum, pretending to be from WIMPs, must be consistent with the 
expectation. However, this shape is exponential, right as it is the case for many background 
sources. 

Obviously, the nuclear-recoil feature, the non-multiple interaction, the uniform event distri- 
bution throughout a detector and the shape of the recoil energy spectrum could not be a clear 
"positive signature" of WIMP interactions. One believes that the following three features of 
WIMP-nuclear interaction can serve as a clear "positive signature" . 

The currently most promising, technically reachable and already used (by the DAMA col- 
laboration) "positive signature" is the annual modulation signature. The WIMP flux and its 
average kinetic energy vary annually due to the combined motions of the Earth and Sun relative 
to the galactic center. The impact WIMP energy increases (decreases) when the Earth velocity 
is added to (subtracted from) the velocity of the Sun. The amplitude of the annual modulation 
depends on many factors — the details of the halo model, mass of the WIMP, the year-averaged 
rate (or total WIMP-nuclear cross sections), etc. In general the expected modulation amplitude 
is rather small (see for example, [26, 27] and [46, 47]) and to observe it one needs huge (at best 
ton scale) detectors which can continuously operate over 5-7 years. Of course, to reliably use 
this signature one should prove the absence of annually-modulated backgrounds. One should, 
however, also be aware that seasonal modulation can also originate from other scenarios such 
as caustic rings of axions or neutralinos in the halo dark matter distribution [163, 164]. 

Another potentially promising positive WIMP signature is connected with the possibility 
to measure the direction of the recoil nuclei induced by a WIMP. In these directional recoil 
experiments one plans to measure the correlation of the event rate with the Sun's motion (see 
for example, [87, 161, 165]). Unfortunately, the task is extremely complicated (see for example, 
[81-85]). 

The third well-known potentially useful positive WIMP signature is connected with the 
coherence of the WIMP-nucleus spin-independent interaction. Due to a rather low momentum 
transfer a WIMP coherently scatters on the whole target nucleus and the elastic cross-section 
of this interaction should be proportional to A 2 , where A is the atomic number of the target 
nucleus. Contrary to the ^-behavior, the cross-section of neutron scattering on nuclei (due 
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to the strong nature of this interaction) is proportional to the geometrical cross-section of the 
target nucleus (A 2 / 3 - dependence). To reliably use this A 2 -signature one has to satisfy at least 
two conditions. First, one should be sure that the spin-independent WIMP-nuclear interaction 
indeed dominates over the relevant spin-dependent interaction. This is far from being obvious 
(see for example, [51, 69, 90, 166, 167]). Second, one should, at least, for two targets with 
different atomic number A rather accurately measure the recoil spectra (in the worst case 
integrated event rates) under the same background conditions. Currently this goal looks far 
from being realizable. 

Developing further the idea of this third signature, one can also consider as a possible extra 
WIMP-signature an observation of the similarity (or coherent behavior) of measured spectra 
at different (also non-zero spin) nuclear targets. This possibility relies on rather accurate spin 
structure functions for the experimentally interesting nuclei (see for example, [70, 99]). 

Also in the case of currently very promising event-by-event active background reduction 
techniques (like in the CDMS and EDELWEISS experiments) one inevitably needs clear positive 
WIMP signature (s). Without these signatures one hardly can convince anyone that the final 
spectrum is saturated only by WIMPs. Furthermore with the help of these extra signatures 
one can define the WIMP mass from the spectrum [168, 169]. 

It is known (see for example the discussion in [39] and earlier partly in [124]) that a proof 
of the observation of a dark matter signal is an extremely complicated problem. As pointed 
out above, on this way an interpretation of measurements in the form of exclusion curves helps 
almost nothing. Of course, an exclusion curve is at least something from nothing observed. It 
allows sensitivity comparison of different experiments and therefore allows to decide who at the 
moment is the best 'excluder'. But, for example, supersymmetric theory is, in general, very 
flexible, it has a lot of parameters, and one hardly believes that an exclusion curve can ever 
impose any decisive constraint on it. The situation is much worse due to the already mentioned 
famous nuclear and astrophysical uncertainties involved in the exclusion curves evaluation [86, 
140, 170-175]. This is why, from our point of view, it is not very decisive (or wise) to use very 
refined data and methods (nuclear, astrophysical, numerical, statistical [176], etc) and spend 
big resources fighting only for the best exclusion curve. This fighting could be only accepted, 
perhaps, in the case when one tries to strongly improve the sensitivity of a small detector having 
future plans to use many copies of it in a huge detector array with a total ton-scale mass. 
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As it already has been stressed in [79, 150], in case of a positive DM signal, e.g. the 
detector HDMS has no means to discriminate the signal from background. With a target mass 
of 200 g only, the statistical accuracy within 2 years of measurements is too low in order to 
see the annual modulation, which is nowadays the only available positive signature of WIMP 
interaction with terrestrial matter. The same is completely true for any other potentially very 
accurate low-target-mass direct dark matter search experiment. To have a chance to see the 
annual modulation signature of WIMP-nuclear interaction and to detect dark matter particles, 
as seems to have been done by DAMA, one preferably needs either a GENIUS-like huge setup 
[177, 178] which was planned to operate up to 1000 kg of HPGe detectors of different enrichment 
of 76 Ge and 73 Ge (in a large volume of ultra-pure liquid nitrogen), or, perhaps, a setup with a bit 
smaller mass, which is able to perform permanent data taking over at least several years under 
extremely low background conditions (like for example, the GENIUS-TF experiment [179, 180], 
or a future enlarged EDELWEISS setup). The performing of such experiment seems, however, 
more difficult than originally expected [181-184]. 

VIII. CONCLUSION 

In this review paper the following main questions have been discussed. 

Why do we want to improve the exclusion curves? The answer usually is: to constrain a 
SUSY-like theory. Unfortunately this is an almost hopeless aim due to the huge flexibility 
of such theories and the inevitable necessity of extra information from other SUSY-sensitive 
observables (for example, from LHC, or Tevatron). Almost all experimental groups presenting 
their exclusion curves try to compare them with some SUSY predictions. It is clear from this 
comparison (see for example, Fig. 7) that there are some domains of the SUSY parameter space, 
which are excluded already now by these exclusion curves. What is remarkable, however, is 
that nobody yet has seriously considered — or used otherwise — these constraints for SUSY. In 
short, at the present and foreseeable level of experimental accuracy, simple fighting for the best 
exclusion curve is almost useless, either for real DM detection, or for substantial restrictions 
for SUSY. 

How far one can improve an exclusion curve ? It is almost a question of taste, when one 
should decide to stop speculations on the improvement of the exclusion curve. Almost always 
one can find something to improve the exclusion curve. 



44 



What one would like to see in the future beyond an exclusion curve? New generations of 
dark matter experiments right from their beginning should aim at detection of dark matter 
particles. This will require development of new setups, which will be able to register positive 
signatures of the dark matter particles interactions with nuclear targets. At least the DAMA 
[44-47] and LIBRA [60, 185] experiments are seen on the way. In order to be convincing, 
an eventual WIMP signal should combine more than one of these positive WIMP signatures 
[158, 159]. 

Why one should try to obtain a real recoil energy spectrum? The spectrum allows one to look 
for the annual modulation effect, the only nowadays available positive dark matter signature, 
which can prove existence of dark matter particle interactions with terrestrial nuclei. There are 
also attempts to determine the WIMP mass on the basis of measured recoil spectra [168, 169]. 
Very accurate off-line investigation of the measured spectrum allows one to single out different 
non-WIMP background sources and to perform controllable background subtractions. 

It seems that, at the level of our present knowledge the dark matter problem could not be 
solved independently from other related problems (proof of SUSY, astrophysical dark matter 
properties, etc). Furthermore, due to the huge complexity (technical, physical, astrophysical, 
necessity for positive signatures, etc) to solve the problem of dark matter one should not be 
afraid, but openly use a reliable model-dependent framework — for example the framework of 
SUSY, where the same LSP neutralino should be seen coherently or lead to effects in all available 
experiments (direct and indirect DM searches, rare decays, high-energy searches at LHC, etc). 
Only if such SUSY framework leads to a specific and decisive positive WIMP signature, this 
could mean a proof of SUSY and simultaneous solution of the dark matter problem. It is on 
the other hand absolutely clear, that SUSY although in contrast to others being prefered, since 
requested by 'higher' particle physics theories, such as Superstrings, is not the only candidate 
for the origin of dark matter, and also other scenarios have to be investigated in a comparably 
thorough way. 
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113/679/0-2(R)). The authors thank Dr. I.V. Krivosheina for the long-term and very fruitful 
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